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EXPERIMENTAL INVESTIGATION OF BASE FLOW FIELD AT 
HIGH ALTITUDE FOR A FOUR-ENGINE CLUSTERED 
NOZZLE CONFIGURATION 


SUMMARY 


An experimental investigation of the base flow field for a four-engine 
clustered nozzle configuration was made over the characteristic range of 
center base pressure to ambient pressure ratio range from near 1. 0 to near 
4. 0. Surveys of the base flow were made with static and pitot pressure probes, 
a constant current anemometer, condenser microphones, and wind vanes. The 
reverse flow was confined to the proximity of the model centerline and appeared 
to be an axisymmetric jet decaying from supersonic to subsonic flow. This re- 
verse jet impinged on the center of the heat shield and spread out, forming a 
wall jet. A slowly circulating n separated n flow region was present in the base 
region formed from the viscous pumping of the nozzle jet plumes, the axisym- 
metric reverse jet, and the base wall jet. The static pressure in the circulat- 
ing region was very nearly ambient. Hot-wire measurements indicated that the 
flow of the reverse jet had a relatively low turbulence level. The reverse jet 
originates from the detached shock portion of the ,T inviscid plume-plume impinge- 
ment line" rather than the oblique shock recompression portion; therefore, the 
shear layer type of analysis for calculating the reverse flow is not applicable. 


INTRODUCTION 


The base heating of rocket powered vehicles from the jets of their own 
engines has been a problem of importance for over a decade [1]. The base 
heating is a result of strong radiation from the use of kerosene as a fuel, after- 
burning of turbine exhaust products, and interaction of plumes from multi- 
nozzle configurations. The interaction produces a reverse flow that impinges 
on the vehicle base and results in convective heat transfer. Techniques intro- 
duced by Korst, et al. , [2] and used by Goethert [3] as well as by Page and 
Dixon [4] have been considered to be the major contributions in attempts to 
analyze base flow, but even these are admitted to be limited to qualitative 
estimates. The complexity of the multi-nozzle configurations and the reverse 
flow phenomena has almost eliminated analytical attempts to solve the problem, 
leaving experimental techniques as the only means of obtaining design informa- 
tion. Most experimental investigations have been limited to obtaining surface 



measurements; consequently, the flow field has not been understood. The 
present experimental study examines the base flow itself, so that perhaps 
reasonable analytical models of the base flow can be constructed. 

The center base pressure variation with ambient pressure (altitude) of 
four-engine configurations has become known as the characteristic curve. A 
typical characteristic curve is shown in Figure i. One would have expected 
that the base pressure would simply have decreased with increasing altitude 
until the vent area between the nozzles became choked, and that further de- 
creases in ambient pressure would not have affected the base pressure. Al- 
though this simple model has been known to be incorrect for some time, it has 
remained as a useful concept. The "reversals" that occur in the characteristic 
curve have not been explained. Parametric investigations have shown how the 
characteristic curve varies with heat shield position, nozzle area ratio, nozzle 
spacing, nozzle lip angle, etc. These investigations are reported in the litera- 
ture and some data obtained by R. Matz and Goethert [5] at Arnold Engineering 
Development Center are reproduced in Figure 2. The "choked" conditions 
appear to occur in each case at a nominal base-to-ambient-pressure ratio of 
four. When no interaction between plumes occurs, the ratio is naturally one. 
The present investigation covers the range from near one to near four along the 
characteristic curve. The majority of the data was taken at the three pressure 
ratios indicated in Figure 1 in relation to the characteristic curve; specifically 
these pressure ratios are 20 x iO -4 , 26 x 10' 4 , and 39 x 10~ 4 . The investigation 
was conducted in unheated air, thereby simplifying instrumentation techniques. 
As one would expect, the major features of the reverse flow are similar for hot 
or cold jets. This view is supported in Reference 6 by comparison of the 
characteristic curves obtained from a model test [7] using heated and unheated 
air and from combustion tests [8, 9] of two other four-engine models. If base 
burning were to occur, however, the characteristic flow pattern might be de- 
stroyed. 

The flow field was surveyed in the planes of symmetry, and surface 
pressure measurements were obtained. No external flow existed over the 
model, and there were no configuration variables. 


APPARATUS 

Model 


Details of the model are shown in Figure 3. The nozzles, which are 
conical with an area ratio of 3. 11 and have an exit diameter of 2. 67 inches, are 
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equally spaced on a 6. 6-inch diameter circle. Their external shape is cylin- 
drical. The heat shield is located 2. 0 inches from the nozzle exit plane, giving 
a vent area between nozzles of approximately 2. 0 by 2. 0 inches. This model is 
a larger scale model of the one used to obtain the characteristic curve labeled 
"Matz” in Figure 1. With the heat shield in the nominal position the model 
characteristic curve is given by the one labeled "1 = 0. 8" in Figure 2 (rather 
than "1^ = 2. 0, " because of the scale change) . 


Test Cell and Ejector 


A schematic of the test cell, model, and ejector is shown in Figure 4. 
The dashed line in the ejector represents a second configuration tested by the 
insertion of a cylindrical sleeve. The model and plenum chamber are movable 
with respect to the ejector so that maximum ejector efficiency can be obtained. 
The performance of the first ejector configuration is shown in Figure 5. The 
exhaust plume did not fill this ejector, making the desired pressure ratio of 
10 x io -4 unattainable and making cell pressure control difficult. The cylin- 
drical sleeve was installed to insure precise cell pressure control; however, 
the lowest pressure ratio attainable with the cylindrical ejector was 17 x 10 -4 . 

The test cell ambient pressure was initially measured by a shielded 
impact tube on the cell wall in the nozzle exit plane. At high nozzle plenum 
pressures, the relatively large amount of reversed flow impinging on the cell 
wall in the vicinity of the cell pressure tube produced erratic pressure varia- 
tions. Although most data were taken at a low plenum pressure where this 
behavior did not occur, the cell pressure measurement was switched to a 
nozzle surface location close to the nozzle exit plane and outside the nozzle 
cluster. This pressure, which differs little from the one on the cell wall, is 
representative of the pressure seen by the plumes. Pressure variations on the 
nozzle surface are small, so that the plume sees the same pressure inside the 
cluster as it does outside, at least in the exit plane, or at the nozzle lip. The 
nozzle pressure taps are shown in Figure 6, a photograph taken through the 
access door. The access door mounted on the side of the cell allows entry to 
the cell for frequent instrument changes. 

A three-dimensional traversing mechanism was mounted inside the cell 
above the model for probe manipulation during a run. Operation was from 
either a console or a hand-held remote control box. Two traversing speeds 
were available, designated fast and slow, corresponding to 0. 1 in./sec and 
0. 01 in. /sec. Position potentiometers provided output to digital voltmeters or 
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recording system to an accuracy of 0. 001 inch. An additional traversing 
mechanism providing axial and rotational traversing is located behind the model 
heat shield for probe manipulation through the heat shield center. 


Facility 


The test was conducted in the intermittent tunnel building at the Von 
Karman Facility (VKF) , of Arnold Engineering Development Center, Tullahoma, 
Tennessee. This facility was chosen so that a relatively large model could be 
tested, thereby minimizing traversing probe disturbances to the base flow. 
Schematics of the VKF facility and cell installation are shown in Figures 7 and 
8. Basically, the facility operation is intermittent with dry unheated air sup- 
plied from VKF's 4000 psia high pressure bottle and exhausted into a 200 000 
cu ft vacuum sphere. 

There are two modes of operation: ( 1) The test cell is evacuated by 
opening valve A ( see Figure 8) and then valve B to bring the model plenum to the 
desired total pressure. Valve C is closed. Valve E remains open with vacuum 
pumps operating on the sphere during the run, and the inbleed valve may or may 
not be open, depending on the desired cell pressure. (2) Valve A is opened to 
evacuate the cell. Valve E is then closed. Valve B is open to bring the model 
plenum to the desired pressure. Valve C is opened to allow the vacuum pumps 
to pump on the cell, and the inbleed valve remains closed. Run time is about 
five to six minutes in the first mode and about three minutes in the second. 
Turn-around time depends on the operating plenum pressure; for the range of 
plenum pressures from 20 psia to 100 psia turn-around times vary from 20 
minutes to 50 minutes. This time was used for instrumentation changes. 


DISCUSSION 
Wind Vane 


The small wind vane shown in Figure 9 was very valuable in determining 
the overall characteristics of the base flow field. The vane was slightly damped 
with grease so that it was responsive to flow directions only in the relatively 
high dynamic pressure regions. Probing was restricted to the geometric planes 
of symmetry where the flow component perpendicular to the planes is zero. For 
purposes of discussion, the coordinate system is defined as shown in Figure 10, 
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and will be used in the presentation of data. The origin of the coordinate 
system is located at the heat shield center. Since the XYZ traversing mech- 
anism is mounted above the model, all traverses were made in the vertical 
plane of symmetry, i. e. , the YZ plane. 

It was discovered that the main portion of the reverse flow was confined 
to the centerline of the model — termed n the reverse jet M — and in close prox- 
imity to the heat shield — termed "the wall jet. u Elsewhere, the wind vane was 
not responsive to flow direction and produced erratic results. The flow direc- 
tion throughout the reverse jet was very nearly perpendicular to the heat shield, 
while the flow direction of the wall jet was naturally parallel to the heat shield. 


Static Pressure 


Since the flow through the vent area is confined primarily to the vicinity 
of the heat shield, probes entering the base region through the low dynamic 
pressure region do not significantly alter the base flow field. This was indicated 
by the fact that static pressures on the heat shield and nozzles were unaffected. 
When the probes were moved into the reverse jet, however, the heat shield 
static pressures near the heat shield center did become somewhat erratic 
because of the probe wake. The heat shield pressure distributions for the three 
pressure ratios are given in Figure 11. 

Details of the pressure probe used to obtain static pressures in the base 
region are shown in Figure 12. It is essentially a thin flat circular disk, having 
its edge ( circumference) beveled on the lower surface and having a static orifice 
located at its center. As long as the disk is coincident with the plane of sym- 
metry, i. e. , the YZ plane, the static pressure profiles can be obtained without 
regard for the flow direction; however, since flow directions were known from 
the wind vane, the tube support was always bent so that it was aligned with 
flow, thereby reducing wake disturbances. 

Data were taken in the reverse jet along Y-traverses in 0. 2-inch incre- 
ments for various distances from the heat shield. The data are shown in Figures 
13, 14, and 15 for the three pressure ratios. In general, the static pressure 
gradients in the base region are small. Static pressure distributions in the 
plume impingement regions are more interesting. Continuous traces were 
made in both the Y direction and the Z direction through the near plume region 
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and recorded on a Variplotter. The data were taken at the high traversing speed 
since low speed traverses gave essentially identical results. Since these pres- 
sures could be affected by probe disturbances because of the flow direction down- 
stream of the impingement, the data are shown without scale in composite form 
in Figure 16. A 3- by 3. 5-inch section of the Y-Z plane of symmetry (see 
Figure 16) is bounded on one side by the model centerline (Y = 0) and on 
another side by the nozzle exit plane ( Z = 2) . Also shown is the impingement 
centerline (see Figure 64c) which is in the plane of symmetry and is 2. 33 inches 
from the model centerline. The static pressures exhibit some degree of sym- 
metry about this centerline. It can be seen that the data repeated itself at 
redundant positions. Superimposed on the figure is the inviscid plume impinge- 
ment line for comparison with the peak pressure points. 


Impact Pressure 


Because the impact tubes are not very sensitive to flow direction and the 
wind vane indicated nearly normal reverse flow, impact pressure measurements 
were made with the tube axis parallel to the model centerline. The details of 
the probe are shown in Figure 17. The impact pressure data, taken in the 
same manner as the static pressure data, are shown in Figures 18', 19, and 20 
for the three pressure ratios. The same probe was used in the wall jet with the 
probe axis parallel to the heat shield (the data are shown in Figures 21, 22, and 
23) . Some additional data were taken with a hypodermic tube in order to get 
closer to the heat shield ( see Figures 22 and 23) . 

To obtain a better perspective of the base flow field, the impact data are 
superimposed upon the physical coordinate system in Figures 24, 25, and 26. 

The arrowheads indicate the orientation of the impact tubes relative to the 
model. Also shown in these figures are the inviscid plume impingement lines, 
which will be discussed more fully in a later section. 

The static pressures were combined with the impact data to obtain the 
Mach number profiles shown in Figures 27, 28, and 29. The profiles, along 
the model centerline, are shown in Figure 30 for various pressure ratios. 

These data were taken in 1/2-inch increments. From the figure it can be seen 
that the Mach number increases as the pressure ratio drops and that the position 
of the peak Mach number moves toward the heat shield as does the sonic point. 
Also, there is a marked increase in the static pressure as the reverse jet ap- 
proaches the heat shield at the lower pressure ratio. Evidently the drop in 
Mach number as the flow approaches the heat shield is caused in part by viscous 
dissipation; however, this does not explain the static pressure rise. This static 
pressure rise is puzzling since it begins while the flow is still supersonic. 
Perhaps the presence of smeared-out shocks or multiple shocks near the heat 
shield would explain the data. At the lower pressure ratio, static and impact 
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data were taken continuously on a slow speed traverse along the model center- 
line and recorded on the plotter. There was no evidence from the plot of any 
discontinuity existing in the base region at that pressure ratio; however, it is 
unlikely that the type instrumentation used here would indicate with any clarity 
the presence of small discontinuities were they to occur. Reference 10 attributed 
the smooth decay from supersonic flow to subsonic flow without shocks to viscous 
decay; however, the short running lengths and the increasingly rapid deceleration 
with decreasing pressure ratio, along with the observed rise in static pressure, 
tend to indicate the presence of smeared shocks. Better flow field instrumenta- 
tion, preferably optical instrumentation, is required to determine the true 
nature of the reverse jet recompression. 


Fluctuating Base Pressure 


A Briiel and Kjaer microphone was mounted flush in the center of the 
heat shield to determine if any resonant base cavity pressure fluctuations were 
present. The output of the microphone, recorded on magnetic tape at 60 inches 
per second in the FM mode, was reduced to the power spectral density function 
plots shown in Figures 31 through 34. Data were taken at two plenum pressures, 
20 psia and 70 psia, and at a pressure ratio of 40 x 10~ 4 . A small spike occurs 
in the power spectrum at about 700 cps. Kistler 701A piezoelectric transducers 
mounted off the heat shield center as accelerometers made it apparent that the 
structural vibrations did not affect the microphone output and that a low energy 
resonance of unknown source is present in the base flow cavity. This resonance 
may originate individually or in combination with the plume impingement, the 
plume shear layer, the wind tunnel noise, or possibly from cell cavity resonance. 


Hot Wire Anemometer 


Hot-wire data were taken in the base flow region ( 1) to obtain a degree 
of redundancy in instrumentation and thereby obtain agreement or disagreement 
with the Mach number obtained from the pressure instrumentation; ( 2) to 
attempt to measure adiabatic wire temperatures to be used with pressure data 
to obtain mean flow measurements, and (3) to establish the base flow turbulence 
level. The theoretical methods and the experimental techniques for obtaining 
these data are discussed along with the results in this section. 
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DETERMINATION OF LOCAL FLOW MACH NUMBER 


Method . The hot-wire system provides the heat transfer from the 
wire to the stream and the wire recovery temperature for zero heat transfer. 
The flow field can be determined from this information if one additional thermo- 
dynamic measurement is made and if the relations for heat loss and recovery 
factor are known functions of the flow properties. For the transition region 
between continuum and free molecular flow, theoretical relations are not 
available. Experimental correlations must be relied upon. Dewey’s correla- 
tion [il] which is the best that is presently available, was used in the data 
analysis presented here, where necessary. Dewey's correlation is for an 
infinite wire so that one must take into account the heat transfer from the wire 
to the support. This is done by using the finite-wire equation which is derived 
as follows: 


Starting with the differential heat balance equation for a finite wire 
having internal heat generation, convective heating to the air, and heat loss to 
supports having equal temperatures, one has 


h?rD dz( T - T ) - -f 

w aw 4 dz 


dT 


w 


R 


' w dz 


dz = i 2 


w 


dz 


( 1 ) 


Replacing R byR [i + o;(T -T)] gives the differential equation for the 
w r wr 

wire temperature distribution. The boundary conditions are: 

1. at z = 0, dT /dz = 0 

w 

2. at z = L/2, T = T 

w s 


where z = 0 corresponds to the center of the wire with the assumption of equal 
support temperatures. 

The solution, assuming constant wire thermal conductivity, is 


T 

w 


T - C 
s 


cosh (rnL/2) 


cosh (mz) + C 


( 2 ) 


where 
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( 3 ) 


i 2 R (1 -cxT )/L + h?rD T 
_ r r aw 

hz-D - i 2 R a/L 
r 


and 


m2 = Ir^lT (h7rD - i2R r a / L) • U) 

w 

The mean wire temperature is found by integrating across the wire , so that 


2 L/2 T - C 

T = — f T dz = C + tanh ( v) , (5) 

w L J w v 

where v = m L/2. The finite wire equation can be used in conjunction with 
experimental data to obtain the two-dimensional values of Nusselt number, Nu, 
and adiabatic wire temperature, T . We simply measure the wire temperature 

at any two different wire currents along with the support temperature, which 
then gives two equations and two unknowns — Nu and T . Redundancy can be 

obtained by making measurements at additional currents, thereby giving a check 
on the applicability of the finite wire equation to the physical problem. With 
Nu and T known, the other flow properties can be obtained by any suitable 

iteration scheme. The one used here is as follows: 

1 . Assume the recovery factor to be 1 . 0: 7 ^ = i. 0 

2. Assume then that the Mach number is 1. 0: = 1. 0 

3. Compute stagnation temperature: T q = T /rj 

4. Use the static pressure from the experimental data and the mass 
flow equation to obtain the Reynolds number: 


puD _ 


PDM f 

m rr 

on o 



( 6 ) 
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5. 


Use M and Rn and Dewey's correlation to obtain 77 


D 


6 . 

7. 

8 . 


Increment and iterate above procedure (steps 2 to 6) until 

Use converged values and Dewey's correlation to obtain Nu^ 

‘ Increment and iterate above procedure ( steps i to 8) until 
Nu^ = Nu 


Experiment . A photograph of a typical hot-wire probe used in the base 
region is shown in Figure 35. Design details are shown schematically in Figure 
36. The hot-wire supports were hypodermic needles epoxied into 1/8-inch O. D. 
stainless steel tubing. The needle separation permitted nominal wire lengths of 
30/1000 inch. At each needle tip, a thermocouple junction was formed from 
iron-constantan lead wire which passed through the hollow needle. All of the 
hot wires that were used for the test were soldered onto the thermocouple 
junctions. Care was taken to insure wire voltage/thermocouple voltage isola- 
tion. The probe lead resistance was measured by placing the probe tip in a 
mercury bath and measuring the resistance at the probe connector with a 
Kelvin bridge. Wire flexure and subsequent failure were reduced by applying 
rubber cement to the wire/support juncture. All wires used were platinum - 
10-percent-rhodium of 0. 0001 inch diameter. 


The voltage across the hot wire, probe, and leads was amplified by a 
Dana Model 3420 dc amplifier before being recorded on an X-Y Plotter. The 
current, selected at the current control panel of the Transmetrics Unit, was 
recorded by hand. The position signal from the survey unit potentiometer was 
fed to a Cohu digital voltmeter and to the X-Y Plotter. Since the Plotter was 
spanned to respond with maximum sensitivity for a selected current, current 
changes would drive the plotter pen off scale so that a dc off-set or a bias 
voltage was required. This off-set was obtained from a precision voltage 
source ( Dial-A -Volt model DAV-46D, General Resistance, Inc.) which allowed 
off-sets with an accuracy of 0. 001 mv. The anemometer and recording system 
are shown schematically in Figure 37. During run periods all the equipment was 
left operating continuously and was calibration checked periodically during a 
run shift. 


The hot-wire probe was mounted above the model to the survey unit in 
a vertical position with the i/2-inch probe tip coinciding with the model center- 
line , and the wire oriented parallel to the X-axis so that for traverses on the 
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vertical plane of symmetry (i. e. , the Y-Z plane at X = 0) , the flow will always 
be perpendicular to the wire. Y-traverses across the reverse jet were made 
from the model centerline ( Y = 0) upward. Z-traverses across the wall jet 
were made from the heat shield ( Z = 0) outward. 

Thermocouples were mounted on the lip of one nozzle, in the test cell, 
and in the model plenum chamber. The temperatures were recorded on the 
VKF standard data system, along with the heat shield and nozzle pressure data. 
The nozzle lip temperature continually decreased during the run approaching 
recovery temperature. Also, the plenum temperature slowly decreased be- 
cause of pressure drop in the VKF high pressure bottle. Two cool-down runs 
were usually necessary to insure a nominal temperature variation during sub- 
sequent runs. Nozzle lip and plenum temperature variation during the run 
affect the reverse jet stagnation temperature and have considerable effect on 
the hot-wire measurements. Attempts to measure current-resistance were 
not made until the flow field became thermally steady as indicated by the wire 
support temperature. 

All R-i data were recorded on Y-traverses in the model exit plane 
( Z = 2) . This one traverse was considered sufficient for indicating agreement 
or disagreement with pressure data, and was repeated on several runs so that 
confidence in the hot-wire results could be ascertained. The resistance and 
temperature of the wire were recorded on the run log before and after each run 
so that changes in wire characteristics could be noted and be used as reference 
conditions for the data reduction. 

An example of the data taken from the plotter is reproduced in Figure 38 
from run 907-6. The three currents used were 0. 392 ma, 3. 0 ma, and 4. 0 ma. 
The fourth trace shown is a repeat of the first current. The probe voltage is 
obtained by adding the off-set voltage to that read from the trace. It is noted 
that the first and fourth trace do repeat fairly well. For the first trace at the 
model centerline, Y = 0, the trace reads 3. 04 mv + 10. 2 mv offset or a total 
of 13. 24 mv, while the fourth trace gives 1. 215 mv + 12. 0 mv offset for a total 
of 13. 215 mv, the difference being 0. 025 mv. This corresponds to 0. 064 ohm 
or about 2 degrees Rankine in wire temperature. As seen from Figure 39, the 
support temperature also varied about 2 degrees. The probe used on this par- 
ticular run had a lead resistance of 0. 785 ohm which must be subtracted from 
that obtained from the plotter to obtain wire resistance. 

Figure 40, an example output from the VKF data system, corresponds 
to the data taken immediately before and after each trace. The variations in 
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cell pressure, plenum pressure, plenum temperature, as well as nozzle lip 
temperature during the run, can be directly observed. The nozzle wall tem- 
perature variation from trace i to trace 4 was also 2 degrees. 

Using the data at Y = 0 and the finite wire equation, curves of Nu vs 
Taw were generated for each wire current. The results are plotted in Figure 

41. It is gratifying that the three curves intersect at a single point, indicating 
that the finite wire equation is an adequate description of the physical problem. 
Using the higher and lower wire currents and the procedure outlined in the 
previous section, the Mach number profile was calculated. The result is com- 
pared with the Mach number obtained from the pressure data in Figure 42 for 
run 907-3. The agreement is considered to be excellent; therefore, the Mach 
number used to reduce the remainder of the hot-wire data will be obtained from 
the pressure data. 

MEAN FLOW MEASUREMENTS 

Method. Having established a reasonable confidence level in the Mach 
number obtained from pitot and static pressure probes we can obtain the flow 
properties with only one additional thermodynamic measurement. The measured 
wire adiabatic temperature ( low wire current) can be obtained rather easily in 
comparison to the previous method. Any time hot-wire measurements are 
made, a sustained effort is required and support by people having a good knowl- 
edge of the electronic equipment is needed. Although making adiabatic wire tem- 
perature measurements is one of the easiest to make, it is still necessary to use 
finite wire corrections, the heat transfer correlations are still needed, and an 
iteration is still necessary. However, only a single iteration is required, as 
compared to the double iteration in the previous technique. The iteration is on 
the local stagnation temperature. From an assumed stagnation temperature, 
the flow properties are obtained and a new stagnation temperature is calculated 
from the finite-wire equation. Since the finite-wire equation is used with only 
very low wire currents, equation ( 5) can be simplified by setting i = 0 with the 
result given by 


T = T + 
w aw 


T - T 
s aw 

v 


tanh ( v ) 


(7) 


where 

v = ( L/D) (Nuk /k ) * 
a w 
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Dividing through by T q and replacing T aw / T 0 b Y the recovery factor n and 
solving for gives 


t'T - T tanh (v) 

o 77 - tanh (.j/)] * (ii) 

This iteration procedure closes rapidly, requiring only two or three iterations. 

Experiment . The data were taken in much the same way as outlined in 
the previous section, except that only one current was used (0. 392 ma) . Y- 
traverses were made at four Z stations: 1. 0, 1. 5, 2. 0, 2. 5. Z-traverses 
were made at six Y stations: i. 0, i. 5, 2. 0, 2. 5, 3. 0, 3. 5. Data were 
obtained for the three pressure ratios except for the wall jet where the pressure 
ratio 20 x lo -4 was not obtained. Examples of raw data plots taken directly 
from the X-Y plotters are reproduced in Figures 43 through 48. The voltage 
on these plots must be corrected as before for lead resistance. The short run- 
time did not allow the flow sufficient time to attain complete thermal equilibrium 
before data recording started. Since the local velocity is relatively insensitive 
to changes in local stagnation temperature, attaining complete thermal equili- 
brium was not a necessity. 

Some of the resulting flow properties are presented in Figures 49 
through 54. The range over which the plenum stagnation temperature, the 
nozzle lip temperature, and the test cell ambient temperature varied during 
the particular run is indicated on the local base flow stagnation temperature 
curves. It is noted that the stagnation temperature does not fall in the expected 
range, i. e. , between the nozzle plenum value and the nozzle lip value. Careful 
reassessment of these data and data reduction techniques did not reveal any 
explanation for the apparent discrepancy. 

TURBULENCE MEASUREMENTS 

Method . Morkovin [12] has established the theoretical basis and opera- 
tional procedures to decompose the hot-wire voltage fluctuations into the various 
thermodynamic and kinematic flow fluctuations. With the set of sensing varia- 
bles, u, p, and T , the voltage fluctuations across the wire from mean voltage 

drop, e, as a function of time, is 


T - T 

Ae = -Ae ( 100 - ) - Ae ( 100 U - U ) + Ae ( 100 ■ ) 

p p u u T T 

o o 

(9) 
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or the voltage fluctuation with respect to mode fluctuations is 


Ae = Ae a + Ae r + Ae £ 
<7 t £ 


( 10 ) 


where <r, r, and £ correspond to the entropy, vorticity, and sound modes. The 

two equations serve to define the sensitivity coefficients, Ae , Ae , Ae. , 

■ P ^ t 

Ae m , Ae , and Ae . 

T <7 T 


These coefficients are functions of the electrical system 


and the thermodynamic and kinematic flow properties. The contribution due 
to the sound mode will be assumed negligible for purposes of data reduction. 
In light of the results obtained, this assumption appears to be valid (see 
Reference 12) . The general expression for these coefficients as derived by 
Morkovin are (with the exception of the sound coefficient) 

Ae = Ae + oiAe^ 
a p T 


Ae = /3Ae_, - Ae 
t T u 

o 


( 11 ) 


where a and /3 are functions of Mach number , and where 


Ae = 


eE 


p 100 


Rn 9Nu w Rn dy 

w Nu 9Rn t 77 9Rn 
wr 


Ae = 


eE 


u 100 


A , Rn 9Nu + i _M_ 9Nu f _ _w_ )_1 M_ dji + Rn 
wjNu 9Rn <7 Nu 9M ( jo; r) 9M + r\ 9Rnl 


a eE 

e T 100 
o 


K + A 


i |k - 1 - n + m ^ 
w ) t t Nu 


9Nu 1 M 


9Nu | 


t t Nu 9Rn 2a Nu 9M f 


(12) 


w 1 M 81 Rn dy 

/ 4. m 

r 2a 7 ] 9M t 7 ) 9Rn 
wr t 


The parameters that affect the wire response can be seen directly. Normally, 
these expressions are reduced to simplified form by assuming some of the 
partial derivatives to be small or negligible; in particular, the term (M/Nu) 
(9Nu/9M) is essentially zero at high Reynolds number and Mach number. For 
the present experiment, no such simplifications occur and the full expressions 
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were used. Using Dewey's empirical heat law [ 11 ] , the various partial deriva- 
tives were calculated (see Figures 55 through 58) . Squaring and averaging of 
equation ( 10) yields 

Ae' 2 = Ae 2 cr’ 2 + Ae 2 t’ 2 + 2Ae Ae o’ r' R , ( 13) 

a t a t err 

where the sound mode has been dropped. By dividing through by Ae 2 and using 
Morkovin’s x-y nomenclature, 

x = Ae /Ae ( 14) 

T cr 

y = Ae'/Ae^ 
the equation becomes 

y 2 = a' 2 + 2 x o-t T > R +x 2 t ' 2 , (15) 

OT 

which is the equation of a hyperbola. By taking data at several overheats and 
using a least mean square hyperbola curve fit, the three unknowns, a 1 , t', and 
R can be obtained. 

<JT 

Experiment. All turbulence data were recorded at a single position in 
the reverse flow — on the model centerline at the nozzle exit plane, i. e. , Y = 0, 
Z = 2. The cell was evacuated, and the current and time constant were meas- 
ured at several overheats. After flow was started and the probe moved into 
position, the current required to obtain a given overheat was measured and the 
time constant calculated by slide rule ( using the fact that the time constant is 
inversely proportional to the square of the current for a given overheat) and 
set in. This time constant was checked using a square wave generator on 
several runs and found to be correct. ( The fact that this could be done 
indicated that the turbulence level was not large. ) The wire voltage signal 
from the Shapiro-Edwards amplifier was fed to a Hewlett-Packard RMS meter 
and its 20-second averaging time output fed to a United Systems Corporation 
Digitec dc voltmeter. The digital output was recorded by hand for the several 
different overheats obtained during a run. The RMS noise was recorded at each 
overheat. The overheating parameter a^ was varied from 0. 0 to 0. 4 in approxi- 
mately 0. 04 increments. The amplifier gain and support temperature were 
recorded. The mean flow parameters were taken from the previously obtained 
mean flow measurements , and sensitivity coefficients were calculated by using 
Figures 55 through 58. The various overheat parameters needed to obtain the 
sensitivity coefficients were obtained in the following manner: The R-i data 
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were curve-fitted by using the theoretical finite-wire equation, and the resulting 

equation was differentiated to obtain A and t from their respective defini- 

w wr 

tions. The zero overheat wire resistance and current exhibited some scatter 

from run to run, caused entirely by variation of tunnel stagnation temperature, 

over which we had no control. The finite circuit impedance factor, e, which 

appears in E, was calculated from the known circuit (with the bridge in the 

circuit) as a function of the wire current. The voltage of the hot wire system 

was increased from 25 to 50 volts to reduce the finite circuit impedance factor. 

The data are shown in Figures 58 through 62, where the RMS voltage has been 

corrected for noise and gain. Using the faired curve through the RMS data, the 

x-y fluctuating diagram was calculated, and the results were plotted in Figure 

63. Superimposed is the least mean square hyperbola curve fit having the 

equation y 2 = 50. 7l55x 2 - 50. 9441x + 53. 3067, and therefore a' = 7. 30%, 

t' = 7. 12%, and R =-0.49. 

or 

Analysis of inviscid Plume Impingement 


Since the test was conducted with a particular geometry and particular 
plenum conditions, it is necessary to interpret the results in the light of 
theoretical analysis to prevent their generalization to situations where they are 
invalid. The effect of the various geometric and fluid dynamic variables upon 
the base flow field is determined by their effect upon the plume interaction 
regions and the shock structure existing in these regions. Mentioned previously 
in conjunction with Figures 24, 25, and 26 were the inviscid plume impingement 
lines, which represent the locus of intersection of solid plumes. If the geom- 
etry and flow conditions are such that an oblique shock can turn the flow into the 
plane of symmetry, then the inviscid plume impingement line represents the 
shock attachment line. This flow situation is represented in Figure 64a, b, 
and c. A reference line is drawn in the plane of symmetry through the impinge- 
ment line at the point where streamline Mj makes an angle with the reference 
line of 0j. The impingement line or the shock origin is swept an angle /3 rela- 
tive to the reference line, which is the intersection line of the plane of symmetry 
and the plane containing streamline M*. The element indicated in Figure 64c 
is enlarged in Figure 65. The initial streamline Mj passes through the shock 
and is turned into the plane of symmetry, becoming streamline M 2 . As viewed 
from A in the figure, the shock plane makes an angle G with the plane of sym- 
metry. Likewise, as viewed from B, the plane containing streamline Mj makes 
an angle y with the plane of symmetry. The plane containing both Mj and M 2 , 
labeled plane 1.02 in Figure 66, is perpendicular to the shock plane. M 2 is 
required to be in this normal plane since that plane contains the resultant 
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normal component of and the resultant tangential component of Mj. There- 
fore, the flow deflection angle 6 and the shock wave angle 9 as related to Mj 
through the Rankine-Hugoniot [13] relations must be measured in this same 
normal plane. M 1} 9, and /3 are known for each impingement point as a result 
of geometry. If G were known, then <5 and 9 could also be obtained from geom- 
etry. G can then be varied through the range from zero to some G , and 

m six 

6 versus 9 can be plotted. The points where this geometric curve of 8 versus 

9 crosses the fluid mechanics solution of <5 versus 9 represent the solution 

points. The geometric curves of 8 versus 9 start on the line 9 = 8 for G = G 

max 

and also end on the line 9 = 8 for G = 0. G is the angle the shock makes 

max 

with the plane of symmetry when M 2 coincides with the impingement line; i. e. , 
the plane, containing Mj and perpendicular to the shock, passes through the im- 
pingement line. The geometric curves for several impingement points are 
plotted in Figure 67 along with the fluid dynamic solution. At the first impinge- 
ment point, there are two solutions corresponding to the weak and strong shocks, 
and, likewise, at impingement point 2, except that the two solutions are at two 
flow deflection angles. At the third impingement point, there are two weak shock 
solutions. At point 4 there is only one solution at the point of tangency. Beyond 

point 4, there are no oblique shock solutions. At point 4, 8 = 8 for the 

max 

component of Mj normal to the impingement line. Figure 68 shows 6, 9 , and G 
plotted as they vary along the impingement line. For a clustered nozzle con- 
figuration where the inviscid impingement line is terminated by the model 
centerline before point 4, one could use a Korst-type calculation to calculate 
the mass reversed along the impingement line and integrate from point 1 to the 
model centerline. However, as pointed out by Chu [14] such a configuration 
would require such a large exit Mach number as to preclude its existence for a 
four-engine cluster. Normally, the shock is detached along a portion of the 
impingement line. For the configuration tested here, the shocks could not be 
attached between the model centerline and the tick marks in Figures 24, 25, 
and 26. Therefore, that portion of the inviscid impingement line has little 
meaning. Evidently, the flow is turned by the shock in a lateral direction where 
it becomes trapped and forms a reverse jet. Outside the cluster, this lateral 
flow moves away from the base region and does not affect its environment un- 
less there is an external flow over the vehicle. The strength of the reverse jet 
is dependent upon the percentage of the impingement line having a detached 
shock. Naturally, the lower the external pressure the greater the percentage; 
however, the strength of the reverse jet is dependent upon other factors as well, 
such as the mass flow along the plume boundary which decreases with decreas- 
ing external pressure. When the external pressure becomes very small, its 
influence upon the base flow dwindles, and the plumes are dependent more upon 
the base pressure environment of their own making. 
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At pressure ratios lower than those tested here, the shock will be de- 
tached along the entire impingement length, in which case the reverse flow from 
the closest impingement point (point i) might or might not affect the heat shield 
environment. In this connection, it is noted that this model was run with a fifth 
nozzle mounted at the center of the cluster, and although the recompression 
shocks were detached all along the impingement lines, the heat shield and nozzle 
static taps indicated ambient pressure. This occurred at a pressure ratio of 
26 x 10 -4 . Even though the total mass reversed was greater than that from the 
four-engine configuration at the same pressure ratio, it did not significantly 
affect the base environment. 

There is yet to appear in the literature solutions for the impingement 
regions, in particular, solutions pertaining to the detached shock and lateral 
flow. It appears that time-dependent finite difference techniques have the best 
chance of solving this problem. The laser-doppler [ 15] technique for obtaining 
flow information in these impingement regions is technically feasible at present, 
and perhaps data pertaining to this problem will be generated within the following 
year. 


CONCLUSIONS 


The flow field in the base region of a four-engine cluster consists of an 
almost axisymmetric turbulent jet impinging upon the heat shield center, form- 
ing a wall jet. This pattern persists throughout the "reversal" portion of the 
characteristic curve. The reverse jet originates from the lateral flow emanat- 
ing from behind the detached shock portion of the plume impingement. The 
viscous flow emanating from the attached shock portion of the plume impinge- 
ment region does not significantly affect the base environment and is pumped 
from the base region by the reverse jet and wall jet. 

The reverse jet at first accelerates to supersonic flow and then decel- 
erates as it approaches the heat shield. The deceleration becomes more rapid 
with increasing altitude as the peak Mach number position and the sonic point 
move toward the heat shield. The wall jet boundary layer grows and accelerates 
to sonic velocity in the vicinity of the vent area between nozzles. 
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Heat transfer and pressure distribution over the heat shield probably 
could be predicted for four-engine configurations by using techniques applicable 
to normal jet impingement on flat plates, provided that the strength of the 
reverse jet were known. In the past, calculations of mass reversed have been 
based upon the recompression of the plume viscous mixing layer rather than on 
the actual detached shock-lateral flow pattern that exists. Hopefully, greater 
emphasis will be placed on the real problem in the future. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, December 12, 1968 
933-32-19-00-52 
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FIGURE 6. PHOTOGRAPH OF MODEL AND DIFFUSER 
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FIGURE 7. SCHEMATIC OF TEST SET-UP AT VKF AREA "E" 



FIGURE 8. VALVE LOCATIONS AND NOMENCLATURE 
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FIGURE 9. SCHEMATIC OF WIND VANE 
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FIGURE li. PRESSURE DISTRIBUTIONS ON HEAT SHIELD 
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FIGURE 15. STATIC PRESSURE IN REVERSE JET 
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FIGURE 16. STATIC PRESSURE FIELD IN PLANE OF IMPINGEMENT 
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FIGURE 17. GEOMETRY OF IMPACT PRESSURE TUBE INLET 
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FIGURE 19. IMPACT PRESSURE IN REVERSE JET 
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FIGURE 22. IMPACT PRESSURE IN WALL JET 














FIGURE 25. REVERSE JET AND WALL JET IMPACT PRESSURE PROFILES 
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FIGURE 32 . SPECTRUM OF CONDENSER MICROPHONE OUTPUT LOCATED ON 

CENTER OF HEAT SHIELD 
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FIGURE 35. PHOTOGRAPHS OF HOT WIRE PROBE 
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FIGURE 37. SCHEMATIC OF ANEMOMETER AND RECORDING SYSTEM 
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FIGURE 39. SUPPORT TEMPERATURE DATA FROM PLOTTER 
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(FIGURE 38) 




FIGURE 40. EXAMPLE OUTPUT FROM THE VKF DATA SYSTEM 
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FIGURE 41. FINITE WIRE EQUATION 
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FIGURE 42. COMPARISON OF REVERSE JET MACH NUMBER AS OBTAINED 
FROM HOT WIRE AND PRESSURE INSTRUMENTATION 
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FIGURE 43. HOT WIRE DATA FROM PLOTTER 
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FIGURE 46. SUPPORT TEMPERATURE DATA FROM PLOTTER 
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FIGURE 53. FLOW FIELD PARAMETERS 
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FIGURE 54. FLOW FIELD PARAMETERS 
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FIGURE 55. PARTIAL DERIVATIVES AS OBTAINED FROM DEWEY'S CORRELATION 
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FIGURE 59. EXPERIMENTAL HOT WIRE DATA OBTAINED FROM AN OVERHEATING TRAVERSE 
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FIGURE 60. EXPERIMENTAL HOT WIRE DATA OBTAINED FROM AN 

OVERHEATING TRAVERSE 
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FIGURE 62. SENSITIVITY COEFFICIENTS OBTAINED FROM DEWEY’S 
EMPIRICAL HEAT LAW FORMULATION AND THEORETICAL FINITE 
WIRE OVERHEAT PARAMETERS 
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FIGURE 64. SHOCK ORIGINATING FROM IMPINGEMENT LINE 
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FIGURE 65. NOMENCLATURE FOR IMPINGEMENT ANGLES 
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FIGURE 66. FLOW DEFLECTION AND SHOCK WAVE ANGLES 
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